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Section  1.  Introduction 


Over  the  last  ten  years  there  has  been  a  spectacular 


increase  on  research  and  applications  of  boundary  element 


techniques.  There  has  been  an  explosion  of  books  as  well 


as  a  series  of  Internet io  •  al  Conferences  [10],  specially 


dedicated  to  boundary  element  methods  (BEM).  The  state  of 


the  art  of  asymptotic  error  estimates  of  the  h-version  for 


BEM  is  described  in  several  detailed  articles  (see  for 


example  [24],  [25]).  There  the  theoretical  framework  for 


both  first-kind  and  second-kind  integral  equations  is  the 


theory  of  pseudodifferential  operators.  As  observed  in 


[21]  one  has  for  strongly  elliptic  pseudodifferential 


operators  convergence  of  any  Galerkin  scheme  with  conforming 


boundary  elements;  also  there  holds  quasioptimality  of  the 


Galerkin  error  in  the  energy  norm. 


Almost  all  work  on  BEM  has  been  performed  with  the 


h-version,  where  the  degree  p  of  the  elements  is  fixed. 


usually  on  low  level,  typically  p  =  0,1,2  and  the  accuracy 


is  achieved  by  properly  refining  the  mesh.  Only  recently 


the  p-version  has  been  introduced  into  the  BEM  [1],  [2], 


[3],  [26].  The  p-version  fixes  the  mesh  and  achieves  the 


.v.w  . 


wmmmm 


accuracy  by  increasing  the  degrees  p  of  the  elements 

uniformly  or  selectively.  In  the  finite  element  method 

(FEM)  the  convergence  of  the  p-version  has  been  thoroughly 

investigated  for  one-  and  two-dimensional  boundary  value 

v 

problems  in  a  series  of  papers  by  Babuska  and  others  (4), 

15],  [6],  [7],  [12].  Meanwhile  convergence  results  have  also 

been  derived  for  the  h-p  version  of  the  finite  element 

method  which  is  a  combination  of  the  standard  h-version  and 

the  p-version  ]4],  [8],  [13],  [14]. 

In  this  paper,  we  prove  the  convergence  of  the  p-version 

for  some  Galerkin  boundary  element  schemes  which  use  first- 

kind  integral  equations.  In  Section  2  we  introduce  the 

function  spaces  and  corresponding  norms  used  later  on.  In 

Section  3.1  we  show  that  the  rate  of  convergence  of  the 

p-version  is  an  optimal  one  in  the  H1//2  and  H  1//2-norms 

1  2 

generalizing  known  results  for  H  and  L  -norms.  In  Sections 

3.3  we  approximate  singular  functions  by  the  p-version  in 
~ 1/2  ~-i/2 

the  H  and  H  -norms  and  we  derive  convergence  rates 
which  are  twice  the  rate  of  the  h-version  with  uniform 
mesh.  In  Section  4  we  apply  the  approximation  results  of 
Section  3  to  the  Galerkin  BEM  for  several  integral  equations 
which  are  strongly  elliptic  pseudodifferential  equations. 

As  examples,  we  consider  the  two-dimensional  screen  Neumann 
and  Dirichlet  problems  in  acoustics  where  sharp  regularity 
results  for  the  solutions  are  available  [22],  [23]. 

Furthermore,  we  give  first-kind  boundary  integral  equations 


governing  the  exterior  Dirichlet  and  Neumann  problems  of 
the  three-dimensional  Helmholtz  equation  and  we  present  the 
convergence  rates  for  the  p-version  of  the  corresponding 
boundary  element  Galerkin  schemes. 


Section  2. 


Notation 


Let  T  be  a  simply  connected,  bounded,  smooth,  closed 
curve  in  IR  and  T  be  a  connected  subset  of  T.  By  C  (F)  , 

0  <  k  <  ^  (k  integer)  ,  we  denote  the  space  of  all  functions 
with  continuous  derivatives  of  order  up  to  k  on  f.  The 
Sobolev  spaces  Hs(f)  are  defined  for  s  >  0  to  be  the  restric 
tions  of  HS  +  ^//^(IR^)  to  f  and  for  s  <  0  by  duality, 

Hs  (f)  =  (H-S (f) )  ' 

0  ~  2  ~ 

with  H  (F)  =  L  (I  )  .  These  spaces  are  used  to  define  the 
corresponding  spaces  of  distributions  on  F,  namely,  for  any 
real  s, 

HS(F)  =  {u  e  HS(f):  supp  u  c  f) 

Hs ( F )  -  { u | r :  u  f  HS(f) }. 

The  above  spaces  are  normed  as  follows.  For  u  defined 
on  F,  let  Pu  denote  any  extension  of  u  on  T  and  u*  denote 
the  zero  extension  of  u  on  f .  Then 


u 


hs  <  r ) 


u*  || 


HS(f) 


(2 


Hull  =  inf{  ||  P,u||  :  tu  f  HS  ( f )  } 

HS(D  IIs  (D 


(2 


Note  that  for  s  >  1/2,  s  /  integer  +  1/2,  HS(F)  is  the  usual 
Hg(I')  space  and  for  -1/2  <  s  <  1/2,  Hs(r)  -  HS(F).  For 
s  "  -1/2,  s  /  integer  +  1/2,  HS(F)  -  (H"S(^)),.  We  will  be 
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particularly  interested  in  the  cases  s  =  1/2  and  s  -  -1/2. 

For  s  =  1/2,  the  space  is  also  denoted  by  H1//2(F), 

1  oo 

with  the  equivalent  norm  (see  (18]) 


=  Hull  1/2  + 

H  7  (D 


M  2, -1/2  , |  2 

(1-x  )  u  j | 

Hu(!) 


(2.3) 


where  for  nplicity,  we  have  assumed  T  =  (~1,+1)  (the  qeneral 
case  can  be  treated  by  affine  maps)  and  x  denotes  the  arc 
length.  In  terms  of  duality,  the  following  relations  hold 

H_1/2(D  =  ( H  1/2  (  r )  )  •  ,  ]V1/2(r)  -  (H 1  /2  ( I  )  )  ’ 


Let  f  be  of  length  2tt,  then  HS(F)  may  be  considered  to 


be  spaces  of  2  r -pr>r  iod  ic  functions.  For  u  (  HS(F)  ,  we  may 
thr>n  write 


u(F)  -  v  a.  cos  j P  +  )  b.  sin  if. 

1  ■'  a  J 


so  that  the  H‘  (F)  norm  may  be  equivalently  defined  by 


(2.4) 


H*(F) 


[  Y  a2 ( 1  +  j  2) S  +  Y  b2(l  +  j2)s] 1/2 
j-0  1  j-1  3 


(2.5) 


For  T  a  smooth  open  arc,  we  will  define  P  (I)  to  be  the 

P 

set  of  all  algebraic  polynomials  cf  degree  less  than  or 
equal  to  p  in  s,  the  arc  length  parameter.  P® ( I )  will  denote 

the  subset  of  polynomials  vanishing  at  the  end  points  of  I. 

N 

Let  us  now  subdivide  F  into  N  pieces,  F  -  u  F.,  such 

i~l  1 

that  F^  is  a  smooth  open  arc  with  end  points  A._j,A^ 

(Aq  -  A^)  .  Then  for  p  ?  n ,  S^fT)  will  denote  the  set  of  all 

functions  u  defined  on  T  such  that  the  restriction  u | y  to  T. 

i 


'"*/  1.’  •  1  .  ' '  l'  •J'" «  "f  ,  ‘  .  •' 
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belongs  to  P  (I\).  Moreover,  we  set  for  p  >  1, 

V  (f)  =  S  (f)  n  C(0)  (f)  . 

We  may  assume  that  F  may  be  partitioned  analogously  and 

define  S  (T),  V  (D  as  above.  Then  S^(D,  V^(D  will  denote 
P  P  P  P 

the  subsets  of  functions  that  vanish  at  the  end  points  of  F. 

~  -1/2  ~ 

Note  that  S  (F)  (S  (F))  is  a  subset  of  H  /  (F) 

P  P 

CH~ 1/2 ( D )  while  V  ( T)  (V  (F))  is  a  subset  of  H1/2(f) 

(H1/2(D)  and  Vp(D  is  a  subset  of  H1/2(F). 

So  far  we  have  dealt  with  the  one-dimensional  case. 

We  will  also  be  interested  in  a  simply  connected,  bounded, 

~  3  C 

smooth,  closed  surface  F  r  1R  The  definitions  of  ir  (F) 

are  analogous  to  the  previous  case.  We  now  assume  that  F 

is  partitioned  into  curvilinear  quadrilaterals  and  triangles, 
-  n 

i  .  e.  ,  F  =  u  f .  .  Let  Q  and  T  be  the  reference  square  and 
i-1  1 

trianqle  respectively,  then  =  F.  (Q)  or  F^(T),  where  F^ 
is  a  smooth  bijective  mapping.  We  assume  that  the  inter¬ 
section  of  any  two  disjoint  TVs  is  either  the  empty  set  or 
a  common  vertex  or  a  common  side. 

By  Pp(T)  we  will  denote  the  set  of  all  polynomials  of 

2 

total  deqree  <  p  on  the  triangle  T.  P  (Q)  wiH  denote  the 
set  of  all  polynomials  of  deqree  <  p  in  each  variable  on  Q. 

We  de f i ne 


S  5  ( T )  -  { \i  |  u  |  (F.(F))  r  P  ( T)  if  Fj  is  a  trianqle  and 

p  F  p 

ujj.  (F.(O)  r  Pp(Q)  if  F  .  is  a  quadrilateral)  (2.6) 


and  V  (F)  --  S  (F)  n  c'  '  ( F)  . 
P  P 


(2.7) 


Section  3 


roximation  Theorems 


In  this  section,  we  will  be  interested  in  obtaining 


estimates  for  the  approximation  of  functions  in  HS(f),  HS(F) 


and  Hs ( T)  by  piecewise  polynomials  belonging  to  the  poly¬ 


nomial  subspaces  introduced  in  the  previous  section 


3.1.  Approximation  of  Functions  in  H* 


We  first  present  some  results  for  the  case  when  u,  the 


function  being  approximated  is  known  to  lie  in  Hs.  These 


will  be  used  by  us  in  the  next  section  for  approximating 


problems  on  closed  curves  and  closed  surfaces. 


In  what  follows,  f  will  denote  either  a  closed  curve 


or  a  closed  surface. 


Then  tern  3.1.  Let  y  =  F  a*  T.  Let  u  e  Hs(y),  s  >  1/2.  Then 


fro*  p  -  1,2,...  the*e  ex  (-5*4  u  r  v  (y)  iuch  that 

P  P 


u-u  II  <  Cp 

P  H  7  ( Y ) 


<  rn’ (s_l/2) 


Hs(y) 


(3.1) 


u'lifir  the  emit  ant  C  (  i  independent  o  |(  u  and  p  but  depends 


on  s  and  the  partition  on  y.  Mo  leave* ,  fio*  u  €  HS(F), 


<  rn-<=-l/2)  .1/2 


h  ( r ) 


log  7  p  ||u|| 


HS(D 


(3.2) 


P’ror^.  The  estimate  (3.1)  follows  by  interpolating  the  approxi 


mat  ion  estimates  for  the  p-version  obtained  in  the  and  H* 


norm  (see  [6]).  In  [7],  an  alternative  proof  (for  closed 


curves)  using  Chebyshev  expansions  is  provided  in  Theorem  3.2. 


Moreover,  (3.2)  is  also  proved  in  this  theorem,  the  procedure 
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being  similar  to  our  proof  of  Theorem  3.3  in  Section  3.2. 


□ 


The  above  theorem  provides  estimates  for  the  error  of  the 

best  approximation  in  the  H1,/2  and  H1^2  norms.  The  next  theorem 

-1/2 

provides  estimates  in  the  H  norm.  It  has  been  proved  in  (9J 

for  y  being  a  closed  curve  and  it  is  included  here  for  completeness 

Thee  rnm  3.2.  Let  Y  =  f  on  T ,  u  f  Hs(y)  /  s  >  0.  Then  oh 

p  =  0,1,2,...  thehe  exists  u  c  V  (y)  inch  that 

P  P 


u-u 


P  H-1/2<y) 


Cp 


- Cs+1/2) 


Hull 


HS(y) 


(3.3) 


ivhehe  C  is  a  constant  -independent  o  fa  u  and  p  but  depends  upon 
s  and  the  gh id  on  y . 


Phooj.  Let  Up  f.  Vp  (  y  )  satisfy 


/ u  to  df;  =  / uii)  d£  for  all  w  r  V  (>)• 
Y  P  Y  P 

Then,  with  e  =  u-Up,  we  have  (see  [13]) 


I  e  II  o 

HU(y) 


«  Cp  5  ||  u  | 


HS(y) 


Now,  for  arbitrary  v  (  H  (y)  ,  we  have  by  (3.4)  , 

°<y! 

TIvTT 


f  evd  F,  f  e  ( v-y)  df  j|e||  „  (i  v-y  f|  „ 

Y  _  Y  .  HU(y)  H  (y) 

"  TTvTl  “  ~ 


1  1 

H  (Y) 


h1(y) 


H1  (y) 


(3.4) 


(3.5) 


5  CF  1  llell  o 

HU(y) 


where  y  r  V  (y)  satisfies 


II v-y  II  o  *  cp-1  || v  ||  . 


H  (y) 


H  (Y) 


This  yields 


II e  II  i  ^  Cp  (S  +  1)  ||u||  (3.6 

(HL(y))’  hs(y) 

Interpolating  (3.5),  (3.6)  and  using  the  fact  that 
H~1/2(y)  =  (H1/2(Y))’  -  (H°  (y)  V (Y) ) l/2 
=  (  (H1 { y ) )  ’  ,  H° (y)  )  1/2 


we  obtain  (3.4).  n 

Remark  3.1.  For  y  =  V ,  we  have  Hk(r)  =  Hk(F)  .  For  y  -  f, 
we  have  ||  •  ||  <  ||  *  ||  ,  .  Hence,  in  either  case 

h  '  (r)  h'  '  (D 

(3.3)  yields 


I|U'UP|I  -1/2 


H 


(Y) 


Remark  3.2.  Since  V  (y)  c 

- 

also  hold  for  some  u  e  S 

P  P 


s  Cp  (s+1/2)  ||u||  (3.7 

hr(y) 

Sp(y),  we  see  that  (3.3)  and  (3.7) 
(Y)  . 


Remark  3.3.  So  far  we  have  assumed  that  f  and  F  are  smooth. 
The  above  theorems  may  also  be  modified  to  the  case  when 
T  and  f  are  only  piecewise  smooth. 

~l/2 

3.2.  H  Approxima tion  of  Singular  Functions 

We  are  interested  here  in  approximating  functions  that 
are  defined  on  the  curve  F  and  have  square  root  singularities 
at  the  end  points.  For  simplicity,  we  consider  a  function 
u  defined  on  I  -  [-1,+lJ  by 


where  n  is  a  C  function  satisfying 


y(x)  =  1,  -1  s  x  <  -1/2 

=  0  ,  1/2  <  x  <  1 

~1  /2 

We  consider  the  approximation  of  u  in  the  H  (I)  norm  by 

funct ions  in  P  ( I ) . 

P 

I<et  I  =  (We  may  consider  I  to  be  a  closed  circle.) 

Let  u  be  transformed  to  the  periodic  function  u  on  I  by  the 
mapping  x  -  cos  K ,  i.e.,  u ( C )  =  u(x).  Then  we  see  that 

u(f.)  =  (1+cos  C)1//2  X  ( cos  l )  -  »'7x(cos  U(cos(f,/2))  (3.9) 

The  following  lemma  is  taken  from  f7). 

Lemma  3.1.  ||u||  .  ||u||  1  „  K  o  'i  a  nil  u  <  Hl^2(I). 

H  /Z(I) 

The  main  theorem  of  this  section  is  the  following. 


Theorem  3.3.  Let  u  bp  defined  by  (3.8).  Then  p  =  1,2,... 

there  exist  6  a  pcf  uncmi  aT  u^  in  P  (I)  iac.h  that 

p  p 


Up(*l)  =  u(+l) 


(3.10) 


liu"unlLl/?  <  CP  1  log1/2  p  (3.11) 

P  H  '  (I) 

Proofi.  We  first  consider  the  image  u  of  u,  which  (being  even) 
may  be  written  as 


u  (U 


5"  a,  cos  k£. 


(3.12) 


De  f  i  ne 


u  u  +  u 

P  P 


(3.13) 


whpre  u  is  defined  by 
P 


1  a,  cos  k r. 
«  k 


(3.14) 


and  where  u  is  a  linear  function  such  that  u  satisfies  the 

P 

condition  (3.10),  i.e.,  such  that 


u(*l)  -  (u-u  )(M)  -  (u-u  )  (cos  J  (  + 1 )  )  (3.  IS) 

P  P 

We  now  estimate  the  coefficients  a^  in  (3.12) .  We  have 

"  ,  77  F 

a  =  C  .,r  u  cos  kf  =  C  f  >  (cos  f)  cos  (4)  cos  k£  dC 

K  0  0 

-  C  !  A  (cos  O  (cos  {  (  2k7+1  )  r )  +  cos(  (2-V~->  O  ]df 

0  /.  i 

Here,  C  may  represent  different  constants.  Integrating  by 
parts  gives, 

ak  =  C  |  [  x  (cos  K)  (sin  (  ( ^— )  K )  '  +  sin  (  ( -y--)  f, )  •  j—y)  ]  J 


2k-r  J0 


+  /  X  '  (cos  Osin  C  [sin  (  (— O  *  +  s  in  (  ( U  2jpryl 

So  that  with  x' (-1)  =  0  further  integration  by  parts  yields 


ak  I  5  cj-y  +  I  *  X  '  (cos  [  cos  (  <k  ~  7)  jirn  "  cos  ^  (k  + 

k  0 

4  cos  ( (k  - §>  O  2^-t  -  co? nk  + 1>  n^hr] 1 } 


3if)  1 

7)U2kTT 


(  1  17 

"  C\  — y  +  I  /  x"(cos  Osin  K 

"  l,  r\ 


sin  (k  -  -j)  £ 

(Tk+T)  ( 2k-l ) 


sin  (k  +  j)  £, 


sin(k  -  2^  sinfk+^H 
T2k-iV[2k-T)  “  T?Y-TT  m  nr 


<Kl} 
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Now  since  y"  and  all  the  sine  functions  are  bounded  indepen¬ 


dent  of  k,  we  obtain 


(3.16) 


now  estimate  ||u-u  ||  .  .  By  Lemma  3.1,  (2.5),  (3.12), 

P  IT  \ 


and  (3.14),  we  have 


P  '  H1/2  (I) 


H  -  (I) 


lu-u  ||2i/2  ,  =  C  7.  av  ( 1  +k 2  )  1/2 

P  H  7  (I)  k  =  p+ 1 


<  C  Z 

p+1 


(l.k2)1/2 


.  C  C 

which  behaves  like  /  — j  dx  =  — .  Hence, 

p+1  x  p 


I!  u“up  il  1/2  s  £ 

P  H  7  (I)  P 

Next,  we  estimate  ||  u  |j  .  .  Since  u  is  linear, 

H  7  (I) 

II  u  II  1/2  *  C  {  |  u  ( +  1 )  |  +  |u  (-1)  |} 

H  7  (I) 

Now  for  any  x,  by  (3.12),  (3.14), 

cxj  on 

|  (u-u  )  (x)  I  <  E  lak  I  5  Z  "T  -  £ 

p  k=p+l  *  k=p+l  k  p 

Using  (3.15),  (3.19),  and  (3.18),  we  see  that 


"  "  1/2 
HA/  (I) 

which  combined  v;ith  (3.17)  yields 


P  H1^2 ( ! ) 


<  C 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


By  (2.3),  we  know  that 


0 


0 


2-1/2  0 

!  U-Up  II  -1/2  *  H  u_up  li  1/2  4  l|(1~x  >  <u-up>  II  0 

p  HL/Nn  P  w'  d)  P  H  U) 


(3.22) 


Hence  we  must  bound  the  second  term.  We  have 

TT  -A 

+  1  i  l/p  P  f 


2-1  02  ^//p  P  "-  -02  -1 
/  (l-xz)  Mu-u'Vdx  =  (  /  +  /  +  r  ) (u-uu) 1 (sin  r)  Ad; 

p  0  1/p  r  -  ±  P 

P  (3.23) 


1 


Now  —  j-A — j  is  bounded  on  (A  ,  ti  -A).  Hence,  using  (3.19) 


1 

TT - 

P  -  -0  2  -1  c  **  -1 

f  (u-u  )  (sin  FJ  d£  <  /  (sin  £)  dr, 

1/P  P  P  1/p 


1 

IT - 

P 


<  log  p 
P 


(3.24) 


Also,  let  [a]  denote  the  integral  part  of  a.  Then  it 
may  be  verified  that  with 


u  :  =  l  a  .  +  T.  a  cos  £, 

3  2,41 


we  obtain  with  (3.13) 

nr> 

0 


u-u  -  T.  a  (cos  2  j  1 )  +  £  a  (cos  ( 2  j +  1  K~cos  £) 

P  i=(Ei2-l  ’  3 


which  satisfies 


(u-Up)  (cos  *(.♦!))  =  0 


as  required  in  (3.10).  Hence, 
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, 


& 


, 

,v.' 


•< 

fry 


(U-U°) 2  <  C{ ( 


j=(^J 


a2 j (cos  2  j  C~  1 ) ) 


+  ( 


I2  .  +  1  (cos  (2j  +  l)  C-cos  FJ  ) 


<  C{  ( 


•  2.m2 
a2 j  sin  jO 


+  ( 


a2j+1  sin(j  +  l)F,  sin  j£)  } 


so  that 


1/P 


:0,2 


-1 


1/P 


f  (u-u  )  (sin  C)  d£  <  C{  /  ( 


0 


0  j-lEjij 


a2jsin2jO  2  (sin  £)  Xd£ 


1/P 

+  /  ( 


£  aP- .  sin  ( j  +  lK  sin  j^)2(sin  E)  ^f;} 

i-lEJil  3 


0  ;=fP±l 


(3.25) 


Now  sin  j  f  <  j  F, ,  so  that  for  any  e  >  0, 


2  . 


sin'jF,  <  sine  j  ^  sin*"  ^  j  E  s  ( j  C )  ^  •  1  5  (jU 


2-e  . 


Hence,  using  (3.16),  the  first  term  on  the  right  side  of 
(3.25)  is  bounded  by 

2f  -  1/p 

/ 

0 


1/p 

/  ( 


l 


±1,2 

0  j,[P±2j  3 


sin  5  5  2 (1-e ) 


/  t2E  1d5 


C  _-2e  _  C 

7^  p  =  p"5 


The  second  term  may  be  similarly  bounded,  as  may  the  term 


~  ~  0  7  - 1 

/  (u-u^)^(sin  E)  a£ 

1  p 


Using  this  with  (3.23),  (3.24)  gives 


15 


2-1/2  .  0,  ,,  C  log 

(1-x  )  (u-u  )  ||  <  - * 

p  HU(I)  p 


1/2. 


P  P 

which  combined  with  (  3 . 2 1 ) - ( 3 . 22 )  yields  (3.11). 


□ 


In  Section  4  we  will  use  Theorems  3.1  and  3.3  to  bound 
the  error  made  when  a  function  that  is  smooth  in  the  interior 
of  F  and  behaves  like  (3.8)  at  the  end  points  is  approxi¬ 
mated  by  functions  in  V  ( r )  . 

P 

_ i/2 

3.3.  H  Approximation  of  Singular  Functions 

In  this  section,  we  consider  the  approximation  of  func¬ 
tions  u  defined  on  I  =  [-1,+1]  of  the  form 

u  ( x )  =  (x+l)~1/2X(x)  (3.26) 


where  x  is  as  before.  We  are  now  interested  in  approximating 
_ 1/2 

u  in  the  H  (I)  norm  by  functions  in  ( I ) .  To  this  end, 
we  first  prove  the  following  lemma. 


Lemma  3.2. 


let.  f  c 


H1/2 (I) 


Then  f 1 


H  1/2(I) 


and 


IL_ 

H 


1/2 


<  C 


(I) 


H1/2 (I) 


(3.27) 


PAoofi-  Let  4’  e  Cq  ( I )  .  Define  4>*  to  be  the  extension  by  0 
of  ij)  to  ER  .  Then  it  may  be  easily  seen  that  i|»*’  -  ij' '  * ,  so 

that 


ll^  II 


H 


■1/2 


(I) 


-  IN 


•  * 


H 


-1/2 


=  IN'* '  II 


H  1/2 (K) 


<  c  IN-*  1 1 


i  /? 


<  c  IN’  1 1 


i  /? 


(3.28) 


since  ip  e  Cq(I).  (The  inequality  (3.28)  can  be  verified 

taking  Fourier  transforms,  for  instance.) 

~l/2 

Now,  let  f  €  H  (I).  We  use  the  following  definition, 
from  (11],  for  the  H  1//2(I)  norm: 


H  1/2(I) 


<j^c”(I) 


<f '  '<i’>  7 

L  ( I ) 


1  1/2 
H '  (I) 


Hence  with 


<f  •  ,r- 


L2  ( I ) 


=  - 


L2(l) 


We  obtain 


<  f  » i|' '  > 


H  1/2 (I) 


L2(I) 


<J'eC0(I) 


t|.eCQ  ( I ) 


H  '  (I) 

H1/2  (I) 


T*7T  7/2~~ 

H1/Z  (I) 


H  1/2  (I) 


self 


h1/2  m 


by  (3.28).  This  proves  the  lemma, 


With  Lemma  3.2,  we  obtain  the  following  analog  to 
Theorem  3.3. 


T hzo^izm  3.4.  Lzt  u  bz  dzfiinzd  by  (3.26).  Then  thzAz  zxi6t& 

a  polynomial  u  in  P  (I)  6uch  that 
P  P 


I!  „  -1,  1/2 

u-u  L_1  S  Cp  log  p. 


(3.29) 


w  =  /  u  dx  =  2 ( x+ 1 )  x(x)  -  2  /  (x+1)  '  x' (x)dx  =  w.  + 

-1  -1  12 

By  Theorem  3.3,  there  exists  v*  f  P^d)  satisfying 

llwl~Vp!l~l/2  ~  Cp~llog1/2p. 

P  H W  (I) 

Also,  since  \{x)  is  smooth,  lies  in  H2_e  for  any  c  >  0. 

2 

Applying  Theorem  3.1,  there  exists  vp  *  Pp(I)  satisfying 

_  ( 2_  ^  __ 1 \ 

H“2-VpliBl/2(I)  5  CP  ^  1°91/2P  Hw2IIh2_C(i) 

<  Cp  1log1/2p. 

1  2 

Taking  v  =  v  +  v  ,  we  have 
P  P  P 


P  h1/2(i) 


„  „  -1,  1/2 
^  Cp  log  p, 


Finally,  using  Lemma  3.2  and  taking  u  =  v ’ ,  we  obtain  (3.29) 
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Section  4.  The  p-Version  for  Boundary  Elements 

Before  we  apply  the  approximation  results  of  Section  3 
to  the  Galerkin  solutions  of  some  integral  equations  of 
the  first  kind,  let  us  recall  some  basic  facts  on  the 
Galerkin  method.  The  key  to  the  error  analysis  of  Galerkin 's 
method  is  the  following  result  by  Hildebrandt  and  Wienholtz 
[15]  (see  also  (11] ,  [21]) . 

Lemma  4.1.  Let  H  be  a  Hilb  eat  s pace  with  dual  H'  ( not 
necessaai  t  ij  identified  with  H)  and  let  A  be  injective  and 
continuant,  flam  H  into  H'  satisfying  a  Gaadinq  inequality . 

Let  u  (  H  denote  the  Solution  of 

Au  =  f  (4.1) 

whelp  f  e  H'  and  let  u €  S.  c  h  denote  the  Solution  o  (j 

N  N 

the  Gale*,  kin  equations 

<Aun,v>  =  <  f ,  v  >  (4.2) 

Fua  the  amoae  let  (jot  any  4>  e  H  theae  exists  <f  f  S.,  with 

N  N 

4’  -  lim  <J>  <  n  H.  Then  fo\  N  taage  enough  the  Galen  kin 

N N 

equations  (4.2)  aae  uniquely  Solvable  and  theae  holds  with 
a  constant  C  independent  of  u,  u^  and  N  the  eaaoa  estimate 

llu_uNll  <  c  inf  UI«-vN||  :  vN  f  SN)  (4.3) 

whe  a  e  ||*|l  denotes  the  no  am  in  H . 

Next  we  list  several  boundary  value  problems  which  can 


be  reduced  to  strongly  elliptic  integral  equations,  i.e., 

the  correspond i ng  integral  operators  satisfy  a  Gardinq 

inequality  in  appropriate  Sobolev  spaces.  Therefore,  due 

to  Lemma  4.1,  the  corresponding  boundary  element  Galerkin 

methods  converqe  and  the  quas i opt ima 1 i ty  (4.3)  holds  leadinq 

together  with  the  approx i ma t ion  results  of  Section  3  to 

error  estimates  for  the  p-version. 

The  Veum ann  ic’reen  p’trbfcm  in  acoustics  describes  the 

scattering  of  a  plane  wave  at  a  hard  obstacle  F.  Herp  F  is 

given  by  an  oriented  open  arc  being  a  finite  piece  of  a 

2 

smooth  curve  in  IR  .  The  orientation  defines  the  normal 
vector  n  pointing  to  the  side  T 2  (see  Fig.  1).  The  opposite 
side  of  T  will  be  denoted  by  F^.  The  scattering  problem 
leads  to  the  problem:  Find  the  p’ leiiu'ie  amplitude.  6<efd 
u  c  ll2oc(^p)  ia*  ii  faying 

(A+k2)  u  =  0  in  -  IR2  \  f 

3u  _  „  3u  _  . 

3n  gl  •  »n  r2  "  g2  X1  Fig>  l* 

Here  k  /  0,  Im  k  >  0  and  g^,g^  €  H  ^^^(F)  are  given  with 
~-l/2 

g  :=  g^  -  g^  f-  H  (T)  .  In  addition,  we  require  the 
Sommerfeld  radiation  condition 

-  iku  -  o(r  ^2)  and  u  =  0(r  ^2)  as  r  -  |x|  •»  ^ .  (4.5 

» 

From  1 2 3 ]  we  know  that  for  Im  k  >  0,  k  /  0  the  problem  (4.4) 
(4,5)  has  no  e igensolut ions  and  furthermore  it  can  be 
reduced  to  a  hyper s i ngu 1  a r  integral  equation  on  F. 


v  s'*755»a5U  ivv  v 
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Theorem  4.1  [23].  Le£  g^,g2  at1c^  ^  be  p/ven  a  4  above.  Then 

thene  hofdi:  (i)  u  *  H ^oc ( Q r )  6olve&  (4.4),  (4.5)  t rf  and 

~1  /2 

on£«/  (he  jump  [ u ]  |  T.  <-  H  (D  Aatiifiiei  the  inteqial 
i1  q  nation 


D[u]  (z)  :=  -21  [u]  (C) 

r 


4>(z,c)dsr  =  f(z),  z  e  F  (4.6) 

z  C 


«■  <  '  >1 


f  ( z ) 


u’/i e  xe 


=  g1  ( z )  +  g2(z)  +  2/(g1(c)-g2(0  )  jpp  <Mz,C)ds  (4.7) 

F  z  ^ 


4>(z,0  s=  Hq1}  (k  |  z-C  |  ) 


(4.8) 


a  n 


d  Hq  ^  ^  ii  the.  Hanked  function  o  &  the  faiMt  kind  and  oxden 


:pic 


(ii)  Thene  exiiti  exactly  one  solution  <j>  e  H1//2(D, 


v  -  [u]  |  orf  (4.6) 


The  proof  of  the  assertion  (ii)  in  [23]  hinges  on  the 
fact  that  D  is  a  strongly  elliptic  pseudodifferential  operator 
of  order  1.  Therefore  there  exists  a  constant  >  0  and 
a  compact  mapping  C^:  H^2(F)  -*■  H  (F)  such  that 


(4.9) 


Re<  ( D+C.  )  ,  ij'>  5  y  1 1  4'  1 1  _  -i  /  o 

1  ( r )  1  h1/z(F) 

ul/2 


for  every  tp  f  H  (F)  .  This  yields  that  D  is  a  Fredholm 

operator  of  index  zero,  and  the  bijectivity  of  D  follows 

therefore  from  its  injectivity  which  is  guaranteed  by  the 

above  assumptions  on  the  wave  number  k.  Note  that  the 

-1/2 

assumptions  on  g^,g2  imply  f  e  H  (F). 


->w 

■  -  ■f-v- 


JiyJWu  Cuth  i*  fc»  An.  ttafl  .  ■»»« 


Using  localization  and  Mellin  transformation  Stephan 
and  Wpndland  derive  in  [23]  the  following  explicit  regularity 
result  for  the  solution  of  (4.6)  near  the  end  points  z  ,z2 
of 

Lemma  4.2  [23].  Fp\  0  <  o  <  1/2  let  q  ?  H'1/2+'(r),  j  -  1,2, 
be  given.  Then  the  i  elution  [u]|r  ?  H1//2(T)  o  ft  the  integral 
equation  (4.6)  hai  the  j (dm 
2 

[u]  |r  =  l  0!iF,i/2Xi  +  vQ  with  vQ  e  H3//2+n(D,  a,  f  ]R. 

(4.10) 

He\e  denote 6  the  Euclidean  diitance  between  z  <■  F  and  the 
end  point  z  ^  o  |(  F.  x^^  *4  a  C^-cuf-odfl  ((unefton  with 
0  <  Xj_  -  1  and  Xj_  =  1  to  z  ^  =  o  fhe  oppoi/.fe 

end  point,  i  =  1,2. 

The  p-vet4<on  Galeakin  method  for  the  hypersingular 
integral  equation  (4.6)  reads:  Find  v^  »:  V^(D  6uch  that 
with  f  <■  H  3//2(F)  given  by  (4.7)  ((p.’t,  off  <t>  <r  ( r )  the'ie 

ho  Idi 


<  Dv  ,  <t>  >  _  =  <  f ,  d>  >  _ 

P  P  L2(F)  P  l2(F) 


(4.11) 


Here  V  (T)  denotes  the  set  of  continuous,  piecewise 


polynomials  of  degree  <  p  which  vanish  at  the  end  points 


of  F  as  introduced  in  Section  2.  Note  Vp(T)  c 


There  holds  the  following  convergence  result  for  the  Galerkin 


scheme  (4.11) . 
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Theorem  4.2.  Let  p  be  sufficiently  lange .  Then  the  Gale>ikin 
equations  (4.11)  aie  uniquely  solvable  and  foa  the  eaa o* 
between  the  exact  solution  [u]|  f  H1/2(r)  of  (4.6)  and  the 
Gal  e  ik i n  solution  v^  f  ( r )  we  have 


(u1_vd||-i/2  -  cp  1  log1/2p 

*  H  '  (T) 


(4.12) 


whole  the  constant  C  is  independent  of  p. 


Pacof.  We  observe  that  the  operator  D  in  (4.6)  fulfills  the 

requirements  on  A  in  Lemma  4.1  with  H  :  =  H1//2(r),  and 
-1/2 

H  =  H  (F)  since  D  satisfies  the  Girding  inequality  (4.9) 
and  D  is  bijective  from  H1//2(T)  onto  H_1//2(D  due  to 
Theorem  4.2.  On  the  other  hand  {V^fD}  is  a  sequence  of 
approximating  subspaces  of  H1//2(D  as  p  •+  »>  and  therefore 
Vp(f)  is  a  candidate  for  the  subspace  in  Lemma  4.1  with 
p  instead  of  N.  Thus  the  convergence  of  the  p-version  for 
the  Galerkin  procedure  (4.11)  is  an  immediate  consequence 
of  Lemma  4.1.  The  rate  of  convergence  in  (4.12)  follows 
from  the  quasioptimality  (4.3)  together  with  the  regularity 
result  (4.10),  where  the  approximation  result  (3.11)  is  used  to 
approximate  the  singular  part  in  (4.10)  and  Theorem  3.1  is  used 
to  approximate  the  regular  part  v^ .  □ 

Remark  4.1.  (i)  The  decomposition  (4.10)  shows  that  the  exact 

solution  \p  =  [u]  |  p  of  the  integral  equation  (4.6)  belongs 
to  H ^  r'(D  for  any  e  >  0.  Therefore  the  h-version  of  the 
Galerkin  procedure  for  equation  (4.6)  gives  only  an  estimate 
of  order  0(h^2)  for  the  Galerkin  error,  if  a  uniform  mesh  is  used, 
(ii)  Application  of  the  estimate  (3.2)  to  the  quasioptimality 


i  1 1 1; i ii . rn u <i ■■■■>«  i .. ... .  . m . .,. ,7. 
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estimate  (4.3)  qives  with  U'  €  e(F) 


II  n  1 1  ~  1  /  2 

p  H  (D 


„  -1/2+r  ,  1/2 

<  Cp  log  p 


H  J  (  F) 


The  better  estimate  (4.10)  follows  from  Theorem  3.3. 


The  Dirichlet  icfie.cn  problem  in  acoustics  describes 
the  scatterinq  of  a  plane  wave  at  a  soft  obstacle  F.  With 
T  beinq  an  open  arc  as  introduced  above  the  scatterinq  prob¬ 
lem  becomes:  find  the  pressure  amplitude  faield  u  t  (51„) 

satisfying 

(A+k^)u  =  0  in  ftp  =  IR^  \  f ,  u  =  g  on  1  (4.13) 

together  with  the  radiation  condition  (4.5)  for  given 
q  f  H ^^(F)  and  k  /  0,  Imk  >  0. 

We  know  from  [22]  that  with  the  above  restrictions  on 
the  wa\re  number  k  the  Dirichlet  problem  (4.13),  (4.5)  has 

no  ei gensolutions .  Furthermore  this  Dirichlet  problem  can 
be  reduced  to  a  weakly  singular  integral  equation  on  F  [22]. 

Th  cole  m  4.3  [22].  Let  g  e  H  ^  ^  ^  {  T )  be  given  and  k/0,  Im  k  >  0 . 
Then  there  holds  :  (i)  u  €  ^loc^F^  halves  (4.13),  (4.5)  if 
and  only  i  f  the  jump  r  H  (T)  satisfies  the  integral 

equati o n 

V[^H](z)  :=  -2/[^](rJt(z,Odsc  =  2g(z),  z  *  F  (4.14) 


where  t  is  given  in  (4.8) 

(ii)  There  exists  exactly  one  solution  </*  e  H  ^^(F),  i {>  -  [ ] 


F 


s  v 


otf  (4.14)  . 

(iii)  let  g  f  H3/2f0(D,  0  •  n  -  1/2,  be 
the  notation  efr  Lemma  4.2  the  Sr  tut  eon  [ 
cl  (4.14)  has  the  j(e  1m 


a  <  v  e  h 

i 


Tn 


H 


Then 

-1/2 


w  i  t  h 

( r ) 


1  In  J 


2 

v 

i-1 


-1/2 

a  .  p  .  v  . 
11  l 


+  fA  weth  i|n  f 


0 


0 


-1/2  +  0  (  f) 


-  1R  (4.15) 


The  proof  of  assertion  (ii)  in  [?2l  uses  that  the 
single  layer  potential  operator  V  is  a  strongly  elliptic 
pseudodi f ferential  operator  of  order  -1.  Therefore  there 

holds  with  a  constant  y2  and  a  compact  operator 

--1/2  ..  1/2  o 

2  '■  ^  (  :  )  v  H  (!')  the  Carding  inequality 


Re MV+C2 )<Ji,*>  2 


l  ( r ) 


h  1/2 ( r  > 


(4.16) 


f nr  any  i  <  H  3 ^2 ( F )  . 


The  decomposition  (4.15)  is  obtained  in  [22 1  by 
localizing  the  weakly  singular  integral  equation  (4.14)  and 
applying  the  Me  1 1 i n  t rans forma t i on .  The  explicit  form 
(4.15)  of  the  solution  near  the  end  points  z^  ,  i  =  1,2, 
allows  us  to  derive  optimal  error  estimates  for  the  Galerkin 
sol ut ion . 

The  p-veii/on  Gafehkin  method  for  the  weakly  singular 
integral  equation  (4.14)  reads:  find  i|’  r  Sp(T)  4uc h  that 
inith  g  ^  H1/2(F) 


VV'p,4'p>  -  "2g,(fp>  aft  <fp  <  SptT)  (4.17) 


Here  Sp(f)  denotes  the  set  of  piecewise  polynomials  of  degree 


<  p  subordinate  to  a  partitioning  of  F  as  introduced  in 

— 1/2 

Sect  ion  2.  Note  S^(r)  r  H  ( I' )  . 

T  (i  e  c  n  e  m  4.4.  Let  p  be  A  u  6  {,  <  c  c  e  n  t  f  u  f  a  n  a  e  .  Then  the  Gat  e  nk  <  n 
e  .)  na  t  <  c  n  ■>  (4.17)  ale  u  n  <  q  u  e  t  u  Ac ( nab  t  e  and  the  ciiei  between 

the  exact  6c  tut  ion  ’*  a  (4.14)  and  the  Gate’ rh/n  Aofution 
ip  f  S  (F)  (4.17)  Aat<A^e^ 

I!  7  ii_  ■,  /->  ‘  Cp  1  loq1/2p  (4.18) 

P  H  A  7  ^  (  r  > 

ie<<n  a  can* fan f  C  i  nde  pc  nden  t  o  p. 

o 

P'icof)-  Due  to  the  Gardinq  inequality  (4.16)  application  of 

— 1/2 

Lemma  4.1  yields  for  the  choices  A  =  V  and  H  -  H  (F)  with 
1/2 

H'  =  H  (F)  the  convergence  of  the  Galerkin  scheme  (4.17)  . 

Note  that  { S ^ ( F )  }  as  introduced  in  Section  2  is  a  sequence 

-1/2 

of  approximating  subspaces  for  (i)  as  p  ♦  00 .  The  esti¬ 

mate  (4.18)  follows  from  the  quasi  opt.  i  ma  1  i  ty  (4.3)  together  with 
the  regularity  result  (4.15)  where  (3.29)  is  used  to  approximate 
the  singular  part  in  (4.15)  and  Theorem  3.2  is  used  to  approxi¬ 
mate  the  regular  part  □ 

The  extend  on  Neumann  (Vinichtct)  pnobfem  in  acoustics 
describes  the  scattering  of  a  plane  wave  at  a  hard  (soft) 
obstacle  P  being  a  bounded  doma  i  n  in  .  For  simplicity 
we  assume  that  the  boundary  F  of  P  is  a  closed,  smooth, 
simply  connected  surface.  Then  the  scattering  problem  leads 
to  the  problem:  Find  u  f  h|oc(IR^\W)  a  a  t  i  A  fiif  i  nq 

2  3  - 

( A+k  )  u  -  0  <H  m  \  <i 


(4.19) 


(4 . 20) 


on 


?n 


g  on  F  ( Neumann ) 


u  =  f  on  r  ( ViAichie  tj  (4.21) 

((o’!  k  f  0,  Im  k  >  0  together  with  the  radiation  condition 

-  iku  =  o(r  ■L),  u  =  0(r  *)  a  6  r  =  [  x  |  ■+  <» .  (4.22) 

Here  we  make  the  general  assumption: 

In  the  extetioA  Neumann  (Vifiichte  t)  pfiobiem  in  1R^\57 
Let  k  be  di^efient  the  eigenvaiuei  o  ^  the 

intetiioK  VLriichlet  (Neumann)  problem  in  9.  (4.23) 


The  restriction  to  the  three-dimensional  case  is  only 
for  simplicity.  Of  course  we  can  derive  analogous  results 
also  for  the  corresponding  2D  problems.  Easy  modi f ica t ions 
of  the  procedure  in  [16],  [20]  lead  directly  to  a  boundary 

integral  equation  method  for  the  Neumann  and  the  Dirichlet 
problem.  One  obtains  immediately  existence  and  uniqueness 
results  analogous  to  Theorems  4.1  and  4.2.  Let  us  first 
consider  again  the  Neumann  problem  (4.19),  (4.20),  (4.22). 


_  1  /  2  ~ 

Theorem  4.5.  Let  g  e  H  (T)  be  given  with  f g  ds  =  0. 

i  r  -j  _ 

Then  the.Ae  hoidi  with  k  ai  above  (i)  u  e  Hi0c^IR  \  Tl)  t>olvei 
(4.19),  (4.20),  (4.22)  ifa  and  onty  if,  u  c  H1//2(F)  &ati6^ie6 
the  integral  equation 


Du  (  2  ) 


-2/uU) 

r 


3n  9n . 


<M  z  ,  £ )  d  s 

7 


g  ( 7. ) , 


z  f  r  (4.24) 


* ( Z , u  := 


(4.25) 


lk  z-c 
e  1 

- 4 71  \  z-r,  | 


q(z)  ~  q(z)  -  2  /  q  (  O  «r  -  -  l(z,Uds 
~  ^  n 


(4.26) 


(ii)  The  ic  cxihth  exact ty  one  ho  tut  ion  u  *  Hi/2(?)  o  ft  (4.24) 

(iii)  Lef  q  *  HS(F),  s  >  -1/2  and  f  be  anabatic.  Then  the 

Q  1  1  ^ 

6o  tut  ion  u  0((  (4.24)  betongh  to  m  (F)  . 


Correspondingly  using  the  direct  approach  of  [20],  [22] 

one  obtains  for  the  Dirichlet  problem  (4.19),  (4.21),  (4.22): 

T  heo-iem  4.6.  Let  f  H  (F)  be  given.  Then  with  k  a*  above 
there  hetdh:  (i)  u  ?  HJocnR3\r2)  hofveh  (4.19),  (4.21),  (4.22 
i  ft  and  ontg  i  i(  f  H  1//2(F)  hatihftieh  the  integiat  equation 

v  ^(Z)  :=•  -2/  ^(C)t(7.,Uds.  -  f  ( z)  ,  z  -r  T  (4.27) 

r 

with  $  ah  in  (4.25)  and 

f(z)  f(z)  +  2 Suit)  t(z,C)dsr  (4.28) 

T 

(ii)  There  exihth  exact ty  one  hot  at  ion  5S  '  H-1/2«n  »t 
(4.27)  . 

(iii)  Let  f  f  Hs(f),  s  >  1/2,  and  f  be  anafyti c.  Then  the 

ho  tut  ion  oft  (4.27)  be  tong h  to  HS  3(f). 

Vroofth  oft  Theoremh  4.5  and  4.6.  For  brevity  we  sketch  only 
the  main  steps.  The  equivalence  (i)  between  the  boundary 


V_V  iA*  JVlvAW  .  a  JV  b.  -  ■  -  *  .Vi.*  jiVW  *  k  ~  -  ,*  »V>V>  *  »Vj*W  •_> *•  *  >  *  »  * 


value  problems  and  the  integral  equations  is  standard  and 
follows  immediately  from  Green's  formula  (see  [16],  [20]). 

The  existence  and  uniqueness  results  (ii)  of  the  solution 
of  the  respective  integral  equation  are  based  on  the  stronq 
ellipticity  of  the  pseudodifferential  operators  D  and  V, 
i.e.,  there  hold  with  constants  Y^*Y2  >  0  and  compact  mappings 
Ci:  H 1  7  2 ( r  )  -  H~1/2(f)  and  C2  :  H “ 1 7  2 ( T )  -  H1/2(f)  the  Garding 

l  nequal i t ies 


Re- (D+C. ) v, v>  _  >  y,  ||v|) 

1  l/(F) 


H1/2(P) 


(4.29) 


Re<(V+C?)4',H»  2  ~  ^  Y?  Ill'll  _1/?  „  (4.30) 

Z  l/(D  *  H  x//(r) 

for  all  v  e  H1/2(f)  and  i}>  £  H-172  (f)  .  Hence 

D:  H1/2(f)  *  H-1/2(f)  and  V:  H_1/2(f)  ->  H1/2(f)  are  Fredholm 

operators  of  index  zero.  Under  the  assumption  (4.23) 
we  have  for  k  /  0,  Im  k  >  0  that  the  integral  equations 
have  no  eigensolutions .  Hence  the  above  mappings  D  and  V 
are  bijective  yielding  assertion  (ii).  The  regularity 
results  (iii)  in  Theorems  4.5  and  4.6  follow  in  a  standard 
way  from  the  ellipticity  of  the  pseudodifferential  operators 
D  and  V  (see  for  example  [19],  [20]).  □ 

Finally,  we  consider  the  Galerkin  equations  for  the 
integral  equations  (4.24)  and  (4.27)  and  show  the  conver¬ 
gence  of  the  p-version. 

The  Gatvfihin  method  (p-veA-iion)  for  the  integral  equa¬ 
tion  (4.24)  reads  with  V  (P)  defined  by  (2.7): 

P 

f / nd  v  f  V  (F)  inch  that  with  g  c  H  ^72(D  given  bu  (4.26) 

p  p  J  J 


the  ne  hold*  all  4* p  e  V  ( f ) 


<  Dv  ,  <t> 


P-  '  P  -2 


=  <9'$, 


L‘(D  P  L2  ( f ) 

Correspondingly,  the  Galen  In n  method  (p-ve  ns i on)  for 

given  by  ( 2 . f ) : 


(4.31) 


the  intearal  eguation  (4.27)  reads  with  ( ( 


Fold  l  e  S  (T)  such  that  auth  f  c  H  ^  ^ 2  ( F )  qu'ftt  bu  (4.28) 
P  P  ' 


then c.  hold 4  (5 o a  £  £  <ip  c  s  ( r ) 


■Vi|>  ,<t» 


_  -  <f  ,< 


p'  P  L2(?)  ’  P  l2(F) 


(4.32) 


Theonem  4.7.  Let  p  be  sufficiently  lange  and  g  f  HS(f),  s  >  -1 


Then  the  Galenkln  equation s  (4.31)  ane  uniquely  solvable. 
.s  +  1 


(D  be  the  exact  solution  of  (4.24)  and  v  e  V  (F) 

P  P 


Let  u  e  H 

be  the  Galenkln  solution  then  we  have  fon  s  >  -1/2 


u-v  ii  „  <  cP-(s+1/2)  iiin  _ 

P  H  7  (D  HS(D 


(4 . 33) 


with  a  constant  C  Independent  of  u,  g  and  p. 


O 

Pnoof.  Obviously,  since  Garding's  inequality  (4.29)  holds 
the  assumptions  of  Lemma  4.1  are  satisfied  if  we  choose 
A  -  D,  H  =  H1/2(f),  H’  =  H~1/2(f)  and  SN  -  V  (f)  c  H1/2(f). 
Note  that  f  is  a  closed,  bounded,  analytic  surface.  Thus 
for  p  large  enough  Lemma  4.1  guarantees  the  unique  solvability 
of  the  Galerkin  equations  (4.31)  and  the  quasioptimal  esti¬ 
mate  for  the  Galerkin  error 


I!  u-v 


H1/2 (f) 


<  C  inf  {  j|  u-w 


H1/2 (f ) 


wp  e  V  (f) }  (4. 34) 


From  Theorem  4.5  (iii)  we  know  that  for  s  >  -1/2  with  a 


constant  C 


u 


H 


s  +  1 


( r ) 


g  II 


HS(?) 


(4.35) 


Therefore  we  can  apply  the  approximation  result  (3.1)  of 
Theorem  3.1  to  (4.34)  and  obtain  (4.33)  by  using  (4.35). 


Thcoiem  4.8.  Let  p  be  Aufficiently  fa' ige  and  f  <r  Hs(f), 


s  >  1/2. 
6o(vabte 

a  ii  d  -j'  f 
P 

s  >  1/2 


Then  the  Gateakin  equation  (4.32)  ate  uniquefu 
Let  €  Hs_1(f)  be  the  exact  lelution  of,  (4.27) 
Sp(D  he  the  Galeahin  Solution,  then  we  have  foi 


II 

11  Th 


Ol 


H  1/2(f) 


s  Cp 


- (s-1/2) 


llfll 


HS(f) 


with  a  constant  C  independent 


f  and  p. 


(4.36) 


Again,  application  of  Lemma  4.1  gives  the  assertion 
if  we  take  A  =  V ,  H  =  H_1/2(f),  H'  -  H1/2(f)  and 
SN  =  S  (?)  r  H  1/2 (F)  since  the  Garding  inequality  (4.30) 
holds.  Note  again  that  V  is  a  closed,  bounded,  analytic 
surface.  From  Theorem  4.6  (iii)  we  know  that  for  s  >  1/2 
with  a  constant  C 


HshII  s-1  ~  5  c  Hf 

Sn  Hs  l(jt} 


HS(f) 


(4.37) 


On  the  other  hand,  Lemma  4.1  yields 


Ih-VU/2, 


<  C  inf 


H  '  (F) 


I  A 

Ith-* 


p  "h  1/2(F) 


<f>  €  Sp(f)  } 


(4.38) 


Therefore  by  applying  Theorem  3.2  and  Remark  3.1  to  (4.38)  we 
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obtain  with  (4.37)  the  desired  estimate  (4.36).  [] 

Remark  4.2.  Theorems  4.7,  4.8  show  that  for  the  p-version, 
the  rate  of  convergence  obtained  depends  only  upon  the  smooth¬ 
ness  of  the  data.  Hence,  when  f  and  q  are  arbitrarily  smooth, 
one  obtains  arbitrarily  high  rates  of  convergence.  This  is 
in  direct  contrast  to  the  h-version,  where  the  rate  of  conver¬ 
gence  depends  in  addition  upon  the  degree  of  polynomials  used 
and  is  therefore  not  very  high  even  for  smooth  solutions. 

Finally,  we  remark  that  results  analogous  to  the  above 
ones  can  be  shown  for  two-dimensional  crack  problems  in 
linear  elasticity,  since  those  problems  can  be  reduced  to 
first  kind  integral  equations  like  (4.6)  or  (4.14)  for  the 
components  of  the  jumps  of  the  displacement  or  traction 
across  the  crack  line  T.  Regularity  results  analogous  to 
(4.10)  and  (4.15)  hold  for  the  solutions  (see  (17],  [22], 

[23]).  Hence  the  corresponding  Galerkin  schemes  (4.11)  and 
(4.17)  will  lead  componentwise  to  error  estimates  like  (4.12) 
and  (4.18)  with  obvious  modifications. 
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